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ABSTRACT 



We have used a coupled time-dependent chemical and dynamical model to investi- 
£T) [ gate the lifetime of the chemical legacy left in the wake of C-type shocks. We concen- 

trate this study on the chemistry of H2O and O2, two molecules which are predicted 
to have abundances that are significantly affected in shock-heated gas. Two models 
are presented: (1) a three-stage model of pre-shock, shocked, and post-shock gas; and 
(2) a Monte-Carlo cloud simulation where we explore the effects of stochastic shock 
activity on molecular gas over a cloud lifetime. For both models we separately exam- 
ine the pure gas-phase chemistry as well as the chemistry including the interactions of 
molecules with grain surfaces. In agreement with previous studies, we find that shock 
velocities in excess of 10 km s _1 are required to convert all of the oxygen not locked 
in CO into H2O before the gas has an opportunity to cool. For pure gas-phase models 
the lifetime of the high water abundances, or "H 2 legacy", in the post-shock gas is 
~ 4 — 7 x 10 5 years, independent of the gas density. A density dependence for the 
lifetime of H2O is found in gas-grain models as the water molecules deplete onto grains 
at the depletion timescale. 

Through the Monte Carlo cloud simulation we demonstrate that the time-average 
abundance of H 2 - the weighted average of the amount of time gas spends in pre- 
shock, shock, and post-shock stages - is a sensitive function of the frequency of shocks. 
Thus we predict that the abundance of H 2 0, and to a lesser extent 2 , can be used 
to trace the history of shock activity in molecular gas. We use previous large-scale 
surveys of molecular outflows to constrain the frequency of 10 km s _1 shocks in regions 
with varying star-formation properties and discuss the observations required to test 
these results. We discuss the post-shock lifetimes for other possible outflow tracers 
(e.g. SiO, CH3OH) and show that the differences between the lifetimes for various 
tracers can produce potentially observable chemical variations between younger and 
older outflows. For gas-grain models we find that the abundance of water-ice on grain 
surfaces can be quite large and is comparable to that observed in molecular clouds. 
This offers a possible alternative method to create water mantles without resorting to 
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grain surface chemistry: gas heating and chemical modification due to a C-type shock 
and subsequent depletion of the gas-phase species onto grain mantles. 

Subject headings: ISM: chemistry; ISM: molecules - stars - shock waves 
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1. Introduction 

The importance of molecular g the mate- 
rial from which stellar and planetary systems are 
made has led numerous investigators to examine 
the chemical processes that combine atoms into 
molecules. Over the past two decades these stud- 
ies have gradually increased in both complex- 
ity and fidelity and can roughly be divided into 
three generations. The first generation chemical 
model solved the chemical rate equations at equi- 
librium, or steady-state, and demonstrated the 
importance of ion-molecule reactions in driving 



the gas-phase chemistry (Herbst & Klemperer 
1973j ). The second generation of models uti- 



lized the later availability of increased comput- 
ing power to study the time evolution of chemical 



abundances ( Prasad fc Huntress 1980 ; Graedel 
Langer, fc Frerking 198^ ). These models, la- 
beled as "pseudo-time dependent" - "pseudo" 
because Hie chemislry evulves with fixed physi- 



callcolidlllohs' - titled obtJei'Ved physical plo-pel'tles 1 

of dense molecular cores (T% ~ ID — 3U K, n# 2 
= 10 4-6 cm -3 ) to show that chemical equilib- 
rium was reached in ~ 10 6-8 years. Pseudo-time 
dependent models represent the most prevalent 
chnmical model and have been quite successful 
in describing the chemistry of quiescent regions 



in both dark and giant molecular cloud cores ( Lee 
t al. 19961; |Bergin et al. 1997|). 



Second generation chemical models operate 
under one simplifying assumption: that the molec- 
ular gas undergoes no dynamical evolution as it 
chemically evolves. However, molecular clouds 
are; certainly dynamically evulving ubjects. The 



tense mass loss, which manifests itself in ener- 
getic winds, bipolar jets, and large-scale outflows 
(c.f. |Lada 1985| ; [Bachiller 19961) . The impact of 



energetic flows on surrounding quiescent gas can 
compress and heat the gas, in some cases pro- 
viding enough energy to overcome endothermic 
barriers of chemical reactions, sputter or destroy 
grains, or even dissociate molecules. The result 
can be a considerably different chemical compo- 
sition in the shocked gas than observed in quies- 
cent material. Since molecules are the primary 
coolants of the gas, any chemical change induced 
by collapse or shocks can also alter the ensuing 
physical evolution of the cloud. To address such 
concerns, a third generation of chemical models 
has been constructed to build on the successes 
of previous generations by combining chemistry 
with dynamics. These models have been prin- 
cipally directed towards investigating the chem- 
istry of core and star formation (c.f 



et al. 19871; IRawlings et al. 1992| ; [Bergin & 





Prasad 


Bcrgin & 



Langer 1997| ) , the physical and chemical struc- 



ture of shocked gas (c.f. |Draine & McKee 19931 ) 
or even complex scenarios in which gas is cycled 
between low and high density through recurrent 
episodes of low-mass star formation ( |Charnley et 



tcr 

'i<f 



al. 1988a| ). 

One goal of coupled astrochemical models is to 
search for and isolate specific molecular species 
that serve as signposts of particular dynami- 
cal events, such as shocks. Coupled models of 
shocked molecular gas have isolated one molecule 
in particular, H2O, which is predicted to form 
in large abundance in shock-heated gas (|Draine 



widespiead uccunence ui supei iheiinai widihs in- 
molecular lines suggests that, overall, the evolu- 
tion of molecular clouds is not simple quiescent 
evolution at a single gas temperature. Moreover, 
the formation of a low- or high-mass star from a 
molecular condensation involves a collapse that 
increases the density by many orders of magni- 
tude. It has also been recognized that the birth 
of a protostar is associated with a period of in- 



Roberge, fc Dalgarno 1983; Kaufmann &: Neufeld 



19964 b); if temperatures in the shocked gas ex- 



ceed 400 K, the endothermic barriers of a few 
key chemical reactions are exceeded and all of 
the available oxygen not locked up in CO will be 
driven into H2O. This led to the assertion that 
water emission is a ubiquitous tracer of shock- 



heated gas QNeufeld fe Melnick 1987] ) . Unfortu- 
nately, due to constraints imposed by the earth's 
atmosphere, the detection of H2O in interstel- 
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lar clouds from ground-based observatories is a 
challenging task. Nevertheless, several studies 
have observed, and detected, transitions of iso- 
topic water (H^O), and even H2O (c.f. Jacq et 



al. 1988 ; Knacke fc Larson 1991 ; ^muidzinas et 



al. 1 1994 ). Quite recently, absorption from warm 
(T gas > 200 K) water vapor has been unambigu- 
ously detected by the Infrared Space Observatory 



(ISO) towards hot stars ( Helmich et al. 1996 ; 
[van Dishoeck &: Hclmich 1996[ ), and in emission 
in HH54 ( Liseau et al. 1996| ). The inferred water 
abundance in these sources is ~ 1 — 6 x 10~ 5 , 
much larger than predicted by chemical mod- 
els run for cold quiescent conditions (x(H20) 
~10 -7 ), and is consistent with water produc- 
tion in warm gas either though high-temperature 
chemistry - as would apply in shock-heated gas 
or in the near vicinity of embedded sources - or 
via evaporation of water-ice mantles. 

An important question yet to be addressed 
by coupled dynamical and chemical models of 
shocked gas in molecular clouds is how long the 
high abundances of water and other molecules 
persist following the passage of a shock. As 
we will demonstrate, the time needed for shock- 
heated gas to return to its pre-shock tempera- 
ture is several orders of magnitude less than the 
time required for the gas to return to its pre- 
shock chemical composition. Thus, the enhanced 
abundances of water and other species could po- 
tentially exist long after the dynamical effects of 
a shock passage are dissipated. Therefore, the 
chemistry inside a cloud is reflective of, and can 
be used to probe, the physical shock history of 
the molecular gas. 

Motivated by these questions, we present the 
results of a coupled dynamical and chemical 
study of molecular gas that is subjected to shocks 
with velocities greater than 10 km s _1 . In partic- 
ular, we follow the time dependence of the chem- 
istry in a shock-heated gas layer as it cools and 
the quiescent time-dependent chemistry is ulti- 
mately re-established. We present models exam- 
ining this evolution using pure gas-phase chem- 



istry as well as chemistry which includes the in- 
teraction of molecules with grain surfaces. 

In Section 2 we discuss the combined chemi- 
cal and dynamical model. In Section 3 we use 
published surveys of molecular outflows to in- 
vestigate the average rate at which shocks with 
a minimally sufficient velocity to affect the wa- 
ter chemistry - 10 km s" 1 - pass a given re- 
gion of a cloud. Section 4 presents the results 
from our combined models along with one exam- 
ple of quiescent chemical evolution. Two models 
are presented: (1) a three-stage model of pre- 
shock, shocked, and post-shock gas; and, (2) a 
Monte-Carlo cloud simulation in which the re- 
sults from Section 3 are used to examine the ef- 
fects of stochastic shock activity on molecular gas 
over a cloud lifetime. Section 5 discusses the im- 
portance of these results on chemical models of 
molecular gas and also reviews the observations 
required to test these models. In Section 6 we 
summarize our results. 

2. Dynamical and Chemical Model 
2.1. Shocks and Chemistry 

Detailed models of magnetohydrodynamic (MHD) 
shocks in dense molecular clouds have been cre- 
ated by a number of authors. For a compre- 
hensive discussion of various models the reader 
is directed to the review in Draine &; McKee 
(1993). Briefly, MHD shocks in interstellar gas 
are divided into two categories: (1) J-type shocks 
where flow variables are discontinuous across the 
shock; and, (2) C-type shocks where the flow 
variables are continuous across the shock front. 
For J-shocks the energy dissipation occurs in 
a thin layer where the temperature increase is 
high enough to destroy molecules. In C-type 
shocks propagating through a magnetic medium 
with low ionization (the case for dense molecu- 
lar gas) the gas is frictionally heated in a thicker 
layer through the process of ambipolar diffusion 
wherein ions tied to the field drift through the 
neutral fluid. The temperatures produced by this 
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Fig. 1. — Maximum attained gas temperature as a function of shock velocity for density of nu 2 = 10 5 
cm~ 3 (taken from Kaufman &z Neufeld 1996b). The slight drop in temperature at ~15 km s _1 is due to 
the onset of H2O formation and cooling. 



process are as high as several thousand Kelvins, 



but the gas remains in molecular form ( Mullan 
1910; praine 198C| ). In addition, the gas com- 
pression in the wake of C-shocks in the region 
of the temperature peak is not exceedingly large 



(c.f. Praine fc Roberge 1982|) ; to first order, the 
neutral density in the gas which will be most af- 
fected by the temperature increase is approxi- 
mately constant. The magnitude of the heating 
depends on the pre-shock density, the strength of 
the magnetic field, and the velocity of the shock. 
For the purposes of this paper we are mostly con- 
cerned with the non-dissociative C-type shocks 
that drive observed molecular outflows. We ac- 
count for the possibility of dissociative shocks in 
the more comprehensive cloud model presented 
in Section 4.4. 

In the models shown here we are solely con- 
cerned with the chemistry in the post-shock gas, 
which is largely driven by the gas temperature, 
molecular hydrogen density, and the visual ex- 
tinction. To simplify the computational task, 
we adopt many of these conditions from detailed 
MHD shock models (e.g., Kaufman and Neufeld 
1996a,b: hereafter KN96a and KN96b), includ- 



ing the variations in these parameters with both 
time and shock velocity. For C-shocks, the post- 
shock H2 density in the region of the tempera- 
ture maximum can be held constant. In addi- 
tion, most shock models assume that the post- 
shock gas is shielded from dissociating radiation, 
an assumption we adopt as well. Thus, the most 
important remaining parameters to be extracted 
from the MHD shock models are the velocity- 
dependent temperature profiles and their time 
evolution. 

KN96b present a self-consistent calculation of 
the heating, cooling, and chemistry of a MHD 
shock wave propagating through dense (n# 2 = 
10 4-6 ' 5 cm -3 ) molecular gas with low (< 10 -7 ) 
ionization. The magnetic field in their model 
is assumed to be perpendicular to the direc- 
tion of shock propagation with a value of B = 
(jih / cm~^) l l 2 fiG. The assumed values for the 
ionization fraction and magnetic field strength 
are consistent with those estimated in molecu- 
lar clouds (e.g. Guelin et al. 1982; Troland & 
Heiles 1986). In Figure 1 we show the maximum 
shock temperature as a function of the shock ve- 
locity (adapted from KN96b). This dependence 
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Fig. 2. — Shock temperature as a function of time for v 6 
M. Kaufman). 



30 km s 1 and n# 2 = 10 5 cm 3 (provided by 



is shown for a pre-shock molecular hydrogen den- 
sity of 10 5 cm -3 , but the maximum post-shock 
temperature is only weakly dependent on the 
density (KN96b). Clearly the gas temperature 
is a strong function of the shock velocity, and is 
essentially a reflection of the strength of the fric- 
tional heating created by ambipolar diffusion. 

Another important parameter required to ex- 
amine the chemistry of post-shock gas is the 



am ount of time the temperature remains at ele- 



vated levels. In Figure 2 we present the temper- 
ature ciS 3j function of time during the passage of 
a 30 km s" 1 shock in a medium with a pre-shock 
density nn 2 = 10 5 cm -3 . These data have been 
kindly provided by M. Kaufman using the code 
described in KN96a,b. The effect of the shock 
is clearly evident as the gas temperature rises 
to greater than 2000 K in less than 100 years. 
After the passage of the shock, the gas rapidly 
cools within a few hundred years. The principal 
coolants of the gas are H2, H2O, and CO. Thus, 
the passage of a 30 km s~ x shock could heat the 
gas from T 



gas 



10-30 K to ~ 2000 K, and this 
gas will return to 30 K within a few hundred 
years. 



2.2. Chemical Model 

In these calculations we have used the chem- 
ical model initially described in Bergin, Langer, 
& Goldsmith (1995) and updated in Bergin &; 
Langer (1997). These papers examined the gas- 
phase chemical evolution including grain sur- 
face molecular depletion and desorption. For 
this calculation we have used both the gas-phase 



UMIST RATE95 chemical network ([Millar et 
"ah 1997| ) and the gas-grain adaption of this 



network described in Bergin & Langer (1997). 
Thus, we separately examine the pure gas-phase 
chemistry and the evolution including grain pro- 
cesses. The gas-grain models include depletion 
onto grain surfaces and desorption via three des- 
orption mechanisms: cosmic-ray-induced heat- 
ing, photodesorption, and simple thermal evapo- 
ration (see Bergin, Langer, fc Goldsmith 1995| ). 



For these models we use the binding energies of 
molecules on an Si02 grain surface. 

For a given calculation the adjustable vari- 
ables are the density of molecular hydrogen, n# 2 , 
the factor by which the external UV field is en- 
hanced above the average galactic interstellar ra- 
diation field, G , the gas temperature, T gas , and 
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the visual extinction to the center of the cloud, 
Ay. For gas-grain models, the dust temperature 
(Tdust) is a l so a variable. The initial chemical 
conditions are taken from Millar et al. (1991), 
and assume that the carbon, silicon, iron, sulfur, 
sodium, potassium, and magnesium are in ion- 
ized form, while the oxygen and nitrogen atoms 
are neutral. All species are assumed to be de- 
pleted relative to their solar abundances, while 
the heavy atoms (e.g., Fe, Mg) are assumed to 
be more severely depleted than carbon, oxygen, 
and nitrogen. 

For chemical models of interstellar shocks, it is 
particularly important that the chemical network 
contains rate coefficients appropriate for high- 
temperature chemistry. To minimize such con- 
cerns, we only discuss in detail results for a few 
simple carbon- and oxygen-bearing species (C, 
CO, HCO+, O, 2 , OH, and H 2 0) whose chem- 
istry is well characterized at both high and low 
temperatures. As mentioned earlier, when tem- 
peratures exceed ~400 K, water is rapidly formed 
in the gas phase through the following series of 
neutral-neutral reactions 



O + H 2 -> OH + H 



OH + H 2 — > H 2 + H 



(Rl) 



(R2) 



QElitzur fe de Jong 19731 ; Pitzur gc Watson 1978| ). 
Both of these reactions possess significant activa- 
tion barriers (E a (Rl) = 3160 K; E a (R2) = 1660 
K). For very high temperatures, molecular oxy- 
ge n will be destroyod by reaction R3 t 



2 + H 2 -> OH + OH 



(R3) 



(E a (R2) = 28190 K), and eventually processed 
into H 2 through reaction R2. We examine in 
detail the temperature dependence and efficiency 
of reactions Rl-3 in the following section. 

The chemical network includes many other 
species beyond the simple species we emphasize 
in this paper and, unless otherwise noted, these 



species have been included in the model runs. 
Because we do not have the same confidence 
in the high-temperature chemical reactions that 
form these species, we do not provide their pre- 
dicted abundances. The UMIST database con- 
tains rate equations for gas temperatures be- 
tween 10 - 50 K (Millar, Farquhar, & Willacy 
1997), a more consistent calculation for other po- 
tentially important shock tracers, such as SiO 
or SO, requires the incorporation of other reac- 
tions such as those estimated for sulphur-bearing 
species in Pineau des Forets et al. (1993). How- 
ever, in Section 5.2 we do discuss the relevance 
these results might have on the abundances and 
post-shock lifetimes of some of these "shock trac- 
ers." 

Careful attention has also been paid to the ac- 
curacy of the low-temperature reaction rates, al- 
though we have sought to minimize this concern 
through the use of the UMIST database. Of par- 
ticular importance for the abundance of water 
at low temperatures are two reactions that have 
been the subject of considerable debate. These 
are the dissociative recombination rate of H^, 
and the branching ratios of the H30 + recombi- 
nation reaction. For H^, the UMIST RATE95 
database has adopted a rate of a(H^) = 3 x 
10 -8 T -o.5 5 from g mith and Spanel (1993). This 

rate has been slightly revised by the recent mea- 
surements of Sundstrom et al. 1994. The adop- 
tion of the more recent value will not alter our 
results. Recently, there have been two studies of 
the H30 + dissociative recombination branching 



ratios ([Williams et al. 1996| ; [Vejby-Christensen 



et al. 1997 ). The branching paths are: 



H 3 0+ + e" 



H + H 2 
OH + H 2 
OH + 2H 
O + H + H 2 



fi (R4) 

f 2 

f 3 

U 



where /1-4 are the branching ratios. Vejby- 
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Christensen et al. (1997) using an ion storage 
ring measured f\ = 0.33, /2 = 0.18, = 0.48, 
and fi = 0.01, while Williams et al. (1996) used 
an ion flow tube to measure f\ = 0.05, fi = 0.36, 
h = 0.29, h = 0.3. The overall OH produc- 
tion predicted by these two studies are in rea- 
sonable agreement, which is important since OH 
reacts with O to form O2. However, the mea- 
surements differ significantly in their production 
of water and atomic oxygen. For the following 
calculations we have used the branching ratios 
estimated by Vejby-Christensen et al. (1997). In 
order to examine the sensitivity of these results 
to this choice, we present one case using the other 
branching ratios. 

3. Shock Timescales in Molecular Clouds 

Before we examine the timescale for the re- 
assertion of quiescent chemistry after the pas- 
sage of a shock, it is useful to examine how of- 
ten molecular material in molecular cores is sub- 
jected to 10 km s _1 shocks. The choice of this 
particular shock velocity is discussed in Section 
4.2. To estimate the rate at which a parcel of 
molecular gas might be subjected to a 10 km s _1 
shock, we have used a method outlined in Mar- 
gulis, Lada, and Snell (1988). Margulis and Lada 
(1986) surveyed a ~ 0.5° x 2° region of the Mon 
OBI molecular cloud in the emission of 12 CO and 
identified nine separate molecular outflows. To 
examine the effect of the molecular outflows on 
the dynamics of the entire molecular core, Mar- 
gulis et al. (1988) estimate the total momentum 
added to the cloud by the outflows. The mo- 
mentum of an outflow is estimated by either: (a) 
integrating directly over the mass x velocity in 
each velocity channel in the line wings, ignoring 
the possible contribution of outflow material hid- 
den by the line core; and/or, (b) assuming that 
the entire outflow is at a constant velocity, which 
is equated to the maximum observed flow veloc- 
ity (c.f. Lada 1985). Using method (a) provides 
a lower limit to the momentum, because it ig- 
nores the inclination of the outflow and, because 



of the difficulty of differentiating the line wings 
from the core, also ignores the mass in the line 
core. In contrast, method (b) is a strict upper 
limit, because it assumes that all of the mass in 
the flow is at the highest observed velocity. 

Using method (a) and the momenta listed in 
Table 2 of Margulis et al. (1988), the total mo- 
mentum generated by the observed outflows in 
Mon OBI is 163 M km s _1 . Momentum con- 
servation requires that 16 M of molecular ma- 
terial will have been swept up when the outflows 
have reached a velocity of 10 km s _1 . Since the 
total mass in the portion of Mon OBI surveyed 
is 3 xlO 4 M (Margulis and Lada 1986), only 
0.05 percent of the total mass will be affected 
by the energetic flows. Thus, in the momentum 
conserving case, it requires ~1900 generations of 
similar star-formation activity for the entire sur- 
veyed area to have experienced a shock of ~10 
km s _1 . Assuming that the average duration of 
the outflow phase is ~ 2 x 10 4 year (Staude and 
Elsasser 1993), then an upper limit to the shock 
timescale, r s , for gas in Mon OBI to be hit by 
10 km s _1 shocks is < 4 x 10 7 years. Because 
the outflow momentum estimates are strict lower 
limits (see discussion of method (a) above), and 
due to the assumption that subsequent genera- 
tions of star-formation activity occurs continu- 
ously, this timescale is an upper limit. Similarly, 
we can use the momenta calculated via method 
(b) to set a lower limit for this timescale. Using 
method (b), Margulis et al. (1988) calculate the 
total momentum to be 883 M Q km s _1 , which is 
3 percent of the total mass of Mon OBI. There- 
fore, ~340 generations of similar star-formation 
activity are required to process the entire core 
by 10 km s _1 shocks, and the lower limit to the 
timescale is r s (Mon OBI) > 7 x 10 6 year. We 
have performed a similar analysis using the ob- 
servations of L1641 in Morgan et al. (1991) and 
find similar results, i.e., r s (L1641)~ 0.6-2.7xl0 7 
years. 

It's important to note that this "cloud-average" 
approach is useful for investigating the poten- 
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Fig. 3. — Time evolution of chemical abundances relative to H2 for quiescent physical conditions. The 
evolution of the H2O abundance is shown for two values of the H^O" 1 " branching ratio (discussed in Section 
3). 



tial of shocks to alter the global composition 
of a cloud. Clearly though, this is an over- 
simplification. In reality, an observer will see por- 
tions of a cloud that are experiencing active star 
formation, for which the local value of t s may 
be <C 10 6 years, while more quiescent regions of 
the same cloud may have local values of t s closer 
to a few times 10 7 yrs, the upper limit quoted 
above. For instance, NGC 1333 is an example of 
a molecular core where a cluster of outflows has 
been detected within a 15' x 15' (1.5 x 1.5 pc) 
area (see Liseau et al. 1998, Warin et al. 1996 
and references therein). Thus, locally the shock 
timescale in NGC 1333 may be below the global 
value of ~ 10 7 years. Another example is the 
Orion Nebular Cluster where several thousand 
stars are estimated to have been formed within 
the past 1-2 million years (Hillenbrand 1997). If 
each new star has an an associated outflow phase, 
it is certainly possible that the shock timescale in 
the dense core associated with this cluster ( OMC- 
1 ) can be reduced to values well below the global 
value. As will be shown, a range in the values of 
t s between < 10 6 and > 10 7 years spans a large 
range in chemical outcomes and may ultimately 



prove to be the cause of strong compositional gra- 
dients within clouds. 

4. Results 

In this section we present the results of the 
coupled dynamical and chemical model of molec- 
ular gas that is subject to shocks with a range of 
velocities. For quiescent conditions, we assume, 
unless otherwise noted, that n# 2 = 10 5 cm -3 , 
Tgas = Tdust = 30 K, G = 1, and Ay = 20 mag. 
These conditions are similar to those observed in 
Giant Molecular Cloud (GMC) Cores, the sites 
of massive star formation ( poldsmith 1987 ). 



4.1. Quiescent Chemistry 

In Figure 3 we present the time dependence of 
the gas-phase chemistry for the quiescent condi- 
tions appropriate in GMC cores. In this figure 
we also show the water abundance calculated us- 
ing two different values for the H30 + dissocia- 
tive recombination branching ratios (see discus- 
sion in Section 2.2). The observed time evolu- 
tion of chemical abundances is quite similar to 



other current gas-phase chemical models (Lee et 
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al. 1996; [Millar et al. 19971 ). Examination of 



the evolution shows that the carbon abundance 
peaks at early times (t ~ 10 4 years for n# 2 = 10 5 
cm -3 ), while for later times all carbon is locked 
into CO. Equilibrium is reached after ~ 3 x 10 6 
years, when the dominant oxygen reservoirs are 
O, CO, and O2. In addition, due to cosmic-ray- 
induced photodissociation of O2 ( [Sternberg fc 
Da ^garno 1995 ), most oxygen remains in atomic 
form, although the total amount of oxygen in O2 
(2 xx(02)) is similar to the abundance of atomic 
oxygen. The steady-state water abundance is dif- 
ferent for the two measured branching ratios; for 
fi = 0.33 the abundance is x(H 2 0) ~ 4.5 x 10~ 7 , 
while for f 1 = 0.05, x(H 2 0) ~ 1.0 x 10~ 7 . For 
all other calculations we use f\ = 0.33. 

4.2. High- Temperature Chemistry 

In Figure 4 we present the dependence of the 
water abundance on both time and temperature, 
while Figure 5 presents a similar plot for O2 • The 
models are for pure gas-phase chemistry where 
we allow the chemical evolution to proceed for 
10 5 years with T gas = 30 K, at which point the 
chemistry is re-initiated and the temperature is 
increased to the values indicated in Figures 4 
and 5. This method allows for appreciable abun- 
dances of molecular species to build up in the gas 
phase at low temperature and makes it easier to 
gauge the effects of high temperature chemistry. 
Thus, t = in these figures refers to time after 
10 5 years of low-temperature or quiescent chem- 
ical evolution when the abundances of water and 
molecular oxygen are x(H20) ~ 4 x 10~ 7 and 
x(02) ~ 3 x 10~ 7 . The results presented in Fig- 
ures 4 and 5 do not depend in detail on the time 
the chemistry is allowed to evolve under quies- 
cent conditions. 

In Figure 4, for T gas > 500 K (correspond- 
ing to shock velocities, v s > 10 km s" 1 ), wa- 
ter accounts for all of the oxygen not bound as 
CO within ~100 years. For temperatures be- 
low ~400 K, the efficiency of reactions Rl and 
R2 beg 111s to rsipidly decline until Tg^g 200 K 



where the high-temperature reactions are unim- 
portant compared to the low-temperature forma- 
tion pathways. In contrast, the abundance of O2 
shown in Figure 5 exhibits the opposite behavior. 
For T gas = 400 K the O2 abundance increases 
above the initial value through normal time- 
dependent chemistry. However, at t ~ 2 x 10 5 
years, the abundance of O2 exhibits a rapid de- 
crease and reaction R3 effectively processes O2 
into H2O. For even higher temperatures the de- 
struction of O2 becomes more efficient, although 
it requires gas temperatures greater than 2000 K, 
corresponding to v s > 26 km s _1 , to destroy all 
O2 within 100 years. 

As demonstrated in Section 2, the post-shock 
gas will not remain at elevated temperatures for 
the entire time shown in Figures 4 and 5. Us- 
ing the cooling information provided in Figure 
2, we see that the time for the post-shock gas 
to cool to below 200 K, what we term the cool- 
ing timescale, is ~ 100 years. Combining the 
information shown in Figures 2 and 4 suggests 
that shocks with velocities > 10 km s _1 will con- 
vert all of the oxygen to water within a cooling 
timescale. The maximum temperatures reached 
behind shocks with velocities less than 10 km 
s _1 will not be high enough to overcome the 
activation barriers and the abundance of water 
should remain relatively unaltered from its qui- 
escent values. This result is similar to that found 
by KN96a,b who found that the maximum shock 
velocity for the conversion of O to H2O is 10 - 
15 km s _1 . Similarly, shock velocities > 26 km 
s _1 are required to destroy all molecular oxygen 
within the cooling timescale. While the exact 
cooling timescale will vary with shock velocity, 
the cooling will be most efficient when the oxy- 
gen is mainly in the form of water - thus, the 
cooling timescale should not vary appreciably as 
long as T gas > 400 K. 

4.3. Three-Stage Chemical Model 

We envisage that the chemical evolution of gas 
in star-forming molecular clouds that have been 
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Fig. 4. — Abundance of H2O as a function of time and gas temperature. The times listed on the abscissa 
are times after ~ 10 5 years of evolution with T gas = 30 K (see Section 4.1). 
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Fig. 5. — Abundance of O2 as a function of time and gas temperature. The times listed on the abscissa 
are times after ~ 10 5 years of evolution with T gas = 30 K (see Section 4.1). 
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Fig. 6. — Schematic demonstrating the physical evolution of the three-stage chemical model of pre-shock, 
shocked, and post-shock gas. In all stages the density, visual extinction, and UV field enhancement are 
kept constant at uh 2 = 10 5 cm~ 3 , Ay = 10 mag, G Q = 1. 



subjected to non-dissociative shocks as one with 
several potentially episodic evolutionary stages. 
The first stage of evolution is simply quiescent 
chemical and physical evolution. The gas starts 
its chemical evolution at low density with mostly 
atomic constituents and proceeds towards higher 
density with an increasing molecular composi- 
tion. At some time the cloud undergoes a period 
of star formation and the protostars alter the lo- 
cal environment through radiation, winds, and 
outflows. The second stage of evolution occurs 
when a parcel of molecular gas is processed by 
a shock. After the passage of a shock, the gas 
rapidly cools, and the pre-shock time-dependent 
chemistry begins to be re-established. Thus the 
post-shock chemistry represents the third evolu- 
tionary stage. For a given parcel of gas, these 
three stages can be repeated, provided that star 
formation is still active in nearby regions (c.f. 
Charnley et al. 1988a,b). 

To account for this picture, we have con- 
structed a three-stage chemical model, shown 
schematically in Figure 6. In Stage (1) we use a 
pseudo-time dependent calculation of the chemi- 



cal evolution in which the gas is initially atomic 
and the physical conditions for quiescent gas are 
those described earlier. This evolution continues 
with constant physical properties until t = 10 5 
years when the gas undergoes a 20 km s _1 non- 
dissociative shock, which heats the gas and raises 
the gas temperature. At the completion of Stage 
(1), Stage (2) is started by re-initializing the 
chemical evolution and raising the gas tempera- 
ture to the maximum value expected for a 20 km 
s _1 shock (T gas w 1000 K; see Figure 1). In this 
fashion, the chemical abundances after 10 5 years 
of quiescent evolution are the initial conditions 
for the second stage. Based on the discussion 
in the previous section, we allow the evolution to 
continue at high temperature for 100 years. After 
100 years, Stage (3) is begun and the chemistry is 
re-initiated with the abundances of 100 years of 
evolution at elevated temperatures providing the 
initial abundances. For the third stage we use the 
original gas temperature (T gas = 30 K) and allow 
the chemistry to re-equilibrate. Through each of 
these stages, pre-shock (Stage 1), shock (Stage 
2), and post-shock (Stage 3), only the temper- 
ature is altered in step-wise fashion, the other 
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physical properties (n# 2 , G Q , Ay) remain con- 
stant. The step-wise time evolution of tempera- 
ture is a simple approximation of the shock tem- 
perature profile shown in Figure 2. As a check 
on these results we have also run a more complex 
multi-stage model using this temperature profile 
and we find that the results which we present in 
the following paragraphs are unaltered. 

4-3.1. Pure Gas-Phase Chemistry 

We first examine the three-stage chemical model 
using only the gas-phase chemical network. Fig- 
ure 7 presents the results from a three-stage 
chemical model of a 20 km s -1 shock impact- 
ing on quiescent gas. In this figure the top panel 
presents the pre-shock chemical evolution (Stage 
1), the middle panel the evolution in shocked gas 
(Stage 2), and the bottom panel the chemistry of 
the post-shock gas layer (Stage 3). Note that the 
temporal range is different for each panel. 

Examining Figure 7 we observe that the chem- 
ical abundances in Stage (1) are quite similar 
to that observed in Figure 3; carbon is slowly 
processed into CO and the abundances H2O and 
O2 steadily increase. In Stage (2) the tempera- 
ture is raised, and the H2O abundance quickly 
rises to incorporate all of the oxygen that is not 
locked into CO. The formation of H2O from O 
is preceded by formation of OH, which accounts 
for the sharp rise and decline in the hydroxyl 
abundance. However, even at 1000 K, the abun- 
dances of O2 and CO are unchanged. In the 
post-shock evolution (Stage 3), the abundances 
of most species do not change for ~ 10 4 years. 
During this time there is a noticeable increase 
in the abundance of atomic oxygen, suggesting 
that even for times < 10 5 years there is some 
destruction of H2O. But, for these early times, 
the abundance of atomic oxygen is more than 
an order of magnitude below that of H2O. Af- 
ter ~ 4 x 10 5 years of low-temperature evolution 
the chemistry dramatically changes as the abun- 
dance of water rapidly drops and the "normal" 
chemistry, with CO, O, and O2 as the dominant 



oxygen reservoirs, is re- asserted. At this time the 
abundance of CO also exhibits a small drop fol- 
lowed by an immediate increase, which is due to 
CO reacting with OH (producing C02). This re- 
action has been measured at low temperatures in 
the laboratory ( Frost, Sharkey, &: Smith 1991| ). 

In Figure 7 the quiescent stage is allowed to 
evolve for 10 5 years, if this stage is allowed to 
evolve for 10 6 years, or even 10 years, the abun- 
dances of a few selected species in the subsequent 
stages would be altered. Examining the full qui- 
escent evolution shown in Figure 3, if the first 
stage evolves beyond 10 5 years, the abundance of 
atomic carbon is lowered, while molecular oxygen 
will be more abundant. If the shock is not strong 
enough to destroy O2, the larger O2 abundance 
would trap oxygen in molecular form and pro- 
duce a lower H2O abundance in the shock and 
post-shock stages. 

For n# 2 = 10 5 cm -3 , the timescale for the 
re-assertion of low-temperature chemistry, r ra , is 
~ 4 — 7 x 10 5 years. In Figure 8 we examine the 
dependence of T ra on density for three-stage mod- 



els run with n H2 = 10 4 , 10 5 , 10 b , 10 7 cm" 



In 



this figure we present only the water abundance 
in the post-shock stage. Even though chemi- 
cal interactions occur with greater frequency at 
higher density, the chemical relaxation timescale 
does not depend on the density. This is a re- 
sult of cosmic-ray driven chemistry. The princi- 
ple mechanism for re-distribution of water back 
to other oxygen-bearing species is through the 
following two reactions: 



H 2 + H+ -> H 3 0+ + H 2 



H2O + HCO" 



H 3 0+ + CO. 



(R5) 



(R6) 



The first channel is important at early times, but 
the second channel dominates at late times. Be- 
cause of the large abundance of water, the abun- 
dances of HCO + and H^~ decrease in the shock 
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Fig. 7. — Evolution of chemical abundances relative to H2 in the pure gas phase three-stage model as 
a function of time for simple carbon- and oxygen-bearing molecules (C, CO, O, OH, O2, H2O, HCO + ). 
Stage 1 (top panel) represents the pre-shock stage with T gas = 30 K, Stage 2 (middle panel) is the shock 
stage with a higher gas temperature (T gas = 1000 K), and Stage 3 (bottom panel) the post-shock stage 
with T gas = 30 K. Note that the time axis for each panel covers a different range and shows the full 
evolution of chemical abundances from one stage to the next. 
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phase. This effect on HCO + is shown in Fig- 
ure 7. However, in the third stage at t > 10 5 
years, the abundance of HCO + and begin 
to increase and water begins to be efficiently 



we have set the dust temperature equal to the 
gas temperature Tdust = Tg as = 30 K. When 
the gas temperature is raised in Stage 2, the 



dust temperature is kept constant (c.f. Draine 



destroyed through the above reactions. Subse- 



Roberge, fc Dalgarno 1983j ). The dominant des- 



quently, the destruction of HCO + and H3 de- 
creases and the process begins to feed itself and 
the water abundance rapidly drops. 



orption mechanism in this model is thermal evap- 
oration from 30 K dust grains. The desorption 
rate is k(H20) evap = 7.2 x 10~ 15 s _1 (see [Bergin 



The, lack of a dependence of T ra on density is 



Langer, fc Goldsmith 1995 ). In this model we 



also an important clue because, while the frac- 
tional abundance of H3" depends on the den- 
sity, the space density of (x(H^)nj^ 2 ) is con- 



stant with density (Lepp, Dalgarno, & Sternberg 
19cQ). Because r ra is tied to the abundance of 
H^, it will depend on the cosmic-ray ionization 
rate. In these calculations we have used an ion- 
ization rate of = 10 17 s -1 ; for an ionization 
rate a factor of 5 higher, the timescale decreases 
approximately by a factor of 5. 

4-3.2. Gas-Phase Chemistry Including Gas-Grain 
Interactions 

Figure 9 presents the three-stage model run 
for n# 2 = 10 5 cm~ 3 , including grain surface molec- 
ular depletion and desorption. For these models 



allow for the depletion of molecular species onto 
grain surfaces and account for some desorption, 
but we do not allow species to react on the 
grain surfaces. For water-ice this assumption is 
not unreasonable because H2O is a stable end- 
point in models of surface chemistry (c.f. Tielens 
and Hagen 1982) and the majority of gas phase 
H2O depletion will occur in the post-shock layer 
when the water abundance has been raised by 
the shock. In addition, icy grain mantles may be 
removed inside shocks (through sputtering) such 
that all of the oxygen will eventually exist in the 
gas phase in the form of water. Thus, water cre- 
ated by the high-temperature shock will deplete 
back onto the mantle directly from the gas phase 
and will not be formed through grain surface re- 
actions. For the three-stage model we do not 
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Fig. 9. — Evolution of chemical abundances in the gas-grain three-stage model as a function of time for 
simple carbon- and oxygen-bearing molecules (C, CO, O, O2, H2O, HCO + ). All abundances are relative 
to H2. This model includes molecular depletion and desorption from grain surfaces and the abundance 
of H2O on the grain surface (H20(gr)) is presented. Stage 1 (top panel) represents the pre-shock stage 
with T gas = 30 K, Stage 2 (middle panel) is the shock stage with a higher gas temperature (T gas = 1000 
K), and Stage 3 (bottom panel) the post-shock stage with T gas = 30 K. Note that the time axis for each 
panel covers a different range and shows the full evolution of chemical abundances from one stage to the 
next. 16 
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Fig. 10. — Evolution of H2O abundance in Stage 3 from the gas-grain model shown for njj 2 = 10 4 , 10 5 , 
and 10 6 cm -3 . Abundances are relative to H2. 



account for any sputtering of the mantle inside 
the shock, but we will account for sputtering in 
the more general model presented later. 

The quiescent evolution shown in the top panel 
of Figure 9 is not appreciably different from that 
seen in the pure gas-phase case. However, there 
is slow but steady accretion of water from the 
gas phase, such that by t = 10 5 years water 
is more abundant on the grain surface than in 
the gas phase. In contrast, due to lower grain 
surface binding energies, and therefore enhanced 
evaporation rates, both O2 and CO remain in 
the gas phase. In Stage 2 the evolution between 
gas-phase and gas-grain models is almost identi- 
cal; all of the oxygen not in CO, or on grains, 
is processed into H2O. In Stage 3 we see the 
largest differences. Here the abundance of gas- 
phase water exhibits a sharp decline at t ~ 10 5 
years, which is less than r ra inferred from pure 
gas-phase chemistry. The decrease in gas-phase 
water abundance is the result of depletion onto 
grains, which has a timescale of ~ 10 5 years for 
this density. As a result, the abundance of wa- 
ter on grains at these times becomes quite large 
~ 10~ 4 . At later times (t > 10 6 years), the abun- 
dance of H2O on the grain surface decreases as 



the water molecules in the grain mantles begin 
to escape via thermal evaporation. Because CO 
and O2 do not deplete, the oxygen chemistry is 
maintained in the gas phase and the overall gas- 
phase water abundance is only a factor of 2 or 
3 below that of the pure gas-phase case. We 
note that similar chemical evolution will be ob- 
served for lower dust temperatures because both 
CO and O2 will not deplete from the gas phase 
(due to cosmic-ray- induced heating), while the 
more tightly bound H2O molecules will deplete 
from the gas phase (see pergin fc Langer 1997 ). 

To examine the dependence of post-shock gas- 
grain chemistry on density, in Figure 10 we 
present Stage 3 for models run at three differ- 
ent densities, n// 2 = 10 , 10 5 ,10 6 cm -3 . This 
figure demonstrates that the inclusion of grains 
places a density dependence on r ra , which is a 
reflection of the depletion timescale for water. 
For nf/ 2 = 10 4 cm -3 , the depletion timescale is 
equal to T ra while for higher densities the deple- 
tion timescale is less than the gas-phase relax- 
ation timescale. Lastly, due to the large abun- 
dance of water trapped on grains, the equilibrium 
chemistry for gas-grain models is never truly re- 
asserted since the chemistry does not reach equi- 
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librium. 

4.4. Monte-Carlo Model of Molecular Cloud 
Evolution 

The results presented in Section 4.3 suggest 
that if the frequency of shock passages is simi- 
lar to r ra , then the water abundance should be 
quite large (~ 10 -4 ). In contrast, if the shock 
timescale is 3> r ra , then the water abundance 
should be close to the lower values expected in 
quiescent gas (~ 10~ 7 ). In this sense, the abun- 
dance of H2O is a measure of the physical shock 
history within a cloud - more so than the dynam- 
ical information contained in line profiles, which 
can be short-lived compared to the lifetime of 
the chemical legacy. Likewise, because O2 is de- 
stroyed in very strong shocks, the abundance of 
O2 may also chart the history of such shocks. 

To examine the utility of using the abundance 
of H2O and O2 to trace the shock history of 
m olecular gas we have created a more general 



spondence between the shock strength (velocity) 
and post-shock temperature is presented in Fig- 
ure 1. If v s > 45 km s _1 we account for the 
possibility of dissociative shocks and destroy all 
molecular abundances. The next cycle is then 
begun with the atomic initial abundances. The 
entire model is continued for a total of ~ 10 8 
years or ~ 10, 000 cycles. 

The computations required to run thousands 
of cycles is too time consuming to use the full 
UMIST database. We therefore constructed a 
smaller network of 26 species and 145 reactions. 
This small network contained all of the important 
high- and low-temperature pathways to form car- 
bon and oxygen species and was tested against a 
sample run that made use of the entire UMIST 
database, with no significant differences found. 

Finally, models of shocked molecular gas have 
found that species can be desorbed from the grain 
mantle, or even the grain core, in shocks through 
sputtering or grain-grain collisions (c.f 



Draine 



ttJHp! nf sWVs in HnnHs using a Mnnt.P-Csrln Roberge, | Dalgarno 198 | Echilk c et al. 1997|; 



simulation. To construct the Monte-Carlo cloud 
model we use the 10 km s _1 shock timescales de- 
termined in Section 3 for massive star-forming 
clouds and the chemical-dynamical model dis- 
cussed in Section 2. In this model we allow the 
chemistry to evolve in steps of 10 4 years. If the 
shock timescale is < t s >= 10 6 years, then dur- 
ing each time interval there is a 1 percent chance 
that the gas will experience a shock. This model 
is similar to the three-stage model, except that 
the time evolution of each stage is 10 4 years. As 
before, when a shock passes (determined by the 
generated random number), it is assumed that 
in all cases the period of elevated temperatures 
is 100 years; for all other times it is assumed 
that T gas = 30 K. As a consequence of the mo- 
mentum driven winds discussed in Section 3, we 
assume that the frequency of shocks faster than 
a given velocity is inversely proportional to that 
velocity. Thus 20 km s _1 shocks occur with a 
frequency that is a factor of two below that for 
shocks faster than 10 km s _1 shocks. The corre- 



Caselli, Hartquist, fc Havnes 1997| ). The amount 
of gas removed depends on the shock velocity and 
differs between various models. Based on the 
results of Caselli et al. (1997) we assume that 
shocks with v s > 15 km s _1 will entirely remove 
the grain mantle. In these models we do not 
explicitly include a sputtering term, but rather 
simulate these effects by raising the dust tem- 
perature to 100 K. This will increase the rate of 
thermal evaporation and rapidly desorb mantle 
species. 

4-4- 1- Pure Gas-Phase Chemistry 

In Figure 11 we show the number of shocks in 
5 km s _1 bins that are generated by the Monte- 
Carlo cloud simulation with < r s >= 10 6 years 
(shaded histogram) and 10 years (filled his- 
togram). From the discussion of shock timescales 
in Section 3, < r s >= 10 6 years is a value appro- 
priate for regions of massive cluster formation 
(e.g., OMC-1), while < t s >= 10 7 years is in- 
ferred for general GMC core material (e.g., Mon 
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Fig. 11. — Number of shocks in 5 km s^ 1 bins from Monte-Carlo model run for ~ 10,000 cycles. Two 
models are shown, t s = 10 6 (shaded histogram) and 10 7 years (filled histogram). 



R2 or L1641). In both cases, the spectrum of 
shocks is heavily weighted towards lower veloci- 
ties, but there are several fast dissociative shocks, 
even for < t s >= 10 7 years. There are a few 
shocks with velocities higher than 80 km s -1 , but 
they are not included in this figure. 

Figure 12 shows the full time dependence of 
the water abundance from the gas-phase Monte- 
Carlo model with < t s >= 10 6 years (top panel) 
and 10 7 years (bottom panel) using the shock 
spectrums shown in the previous figure. For 
< t s >= 10 6 years, the water abundance in- 
creases for the first 10 5 years and at t ~ 10 6 years 
sharply rises as the gas is shock-heated. The 
abundance decreases in the post-shock gas and 
goes through a series of successive shock heating 
and cooling events. The vertical lines that extend 
below x(H20) = 10~ 12 are the result of dissocia- 
tive shocks. Overall, the water abundance varies 
between x(H 2 0) ~ 10~ 6 and 10~ 4 , with an av- 
erage abundance of ~ 10~ 5 . For a longer shock 
timescale (bottom panel), the water abundances 
show similar variations, but the effects are not as 
dramatic. 

We have computed the time-averaged H2O 



and O2 abundances in the model presented in 
Figure 12 and these results are shown in Figure 
13. To create this plot we compute a running av- 
erage as a function of time, weighting the abun- 
dance by the length of time it remains at a given 
value. In this figure we also present the time- 
averages for another Monte-Carlo model run with 
< t s >= 5 x 10 6 years. For < t s >= 10 6 years, 
the time-averaged water abundance converges to 
~ 10 -5 . For other, longer shock timescales, the 
water abundance takes longer to converge to a 
single value, but the overall amount of water is 
clearly dependent on the shock timescale. These 
effects can be quite large for < r s > greater than 
10 years; the average water abundance will be at 
least a factor of five below the value for < t s >= 
10 6 years. It is also notable that even when 
the shock timescale is 10 7 years, the frequency 
of shocks is high enough that the time-averaged 
water abundance is greater than that expected 
for steady-state chemistry in quiescent gas (com- 
pare with Figure 1 where x(H20) ~ 6 x 10~ 7 ). 
In contrast, because the frequency of shocks with 
velocities high enough to destroy O2 within the 
cooling time is low, the abundance of O2 has a 
negligible dependence on the shock timescale. 
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Monte-Carlo Cloud Model 




Fig. 12. — Abundance of H2O relative to H2 shown as a function of time for the Monte-Carlo cloud 
model. The top panel presents the chemical evolution for r s = 10 6 years and the bottom panel for 
t s = 10 7 years. The vertical lines extending to abundances below 1CT 12 are dissociative shocks. 
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Fig. 13. — Time average abundances of H2O and O2 from Monte-Carlo cloud models run with variable 
shock timescales. This model is for pure gas-phase chemistry at nn 2 = 10 5 cm~ 3 and all abundances are 
calculated as relative to H2. The time-averaged abundances are calculated as a running average, weighted 
by the amount of time the abundance remains at a given amount. 
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Fig. 14. — Time averaged H2O and O2 abundances (relative to H2) from Monte-Carlo models run with 
two separate shock timescales and three different densities. 
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Fig. 15. — Same as in Figure 13, but for gas-grain Monte-Carlo model. 



Figure 14 presents the time-averaged abun- 
dances for models run at < t s > = 10 6 and 10 7 
years for three different densities. The results 
shown in this figure confirm the earlier finding 
that for pure gas-phase chemistry and a con- 
stant cosmic-ray ionization rate, the chemical re- 
laxation timescale has little dependence on the 
density. Similarly, the time-averaged H2O abun- 
dance is not dependent on density. While the 
relaxation timescale does not depend on the den- 
sity it does have a dependence on the cosmic-ray 
ionization rate (see discussion in Section 4.3.1). 
The relaxation timescale is inversely proportional 
to the ionization rate, such that a Monte Carlo 
model with a shock timescale t s and cosmic 
ray ionization rate of Qh 2 = 10 17 s_1 would 
roughly apply for a model with a shock timescale 
< t s >'=< t s > /(Ck/lO- 17 *- 1 ), where C'h 2 is 
a different cosmic ray ionization rate. 

Gas-Phase Chemistry Including Gas-Grain 
Interactions 

The time-averaged abundances calculated in 
the Monte-Carlo model, which include interac- 
tions between the gas phase with grain surfaces, 
are shown in Figure 15. These results are pre- 



sented for three densities, n# 2 = 10 4 , 10 5 , 10 6 
cm~ 3 and two shock timescales < t s >= 10 6 and 
10 7 years. As before, species are fully evaporated 
from the grain mantles if shock velocities exceed 
15 km s _1 . The inclusion of depletion onto grains 
introduces a density dependence into the overall 
water abundance. Although the H2O abundance 
does converge to a single time-averaged value, 
this value does have a dependence on the shock 
timescale. Therefore, to compare these results 
with observations, prior information about the 
density is required in order to differentiate be- 
tween the various solutions. For molecular oxy- 
gen, the inclusion of gas-grain interactions intro- 
duces a dependence on the average abundance 
with density (because of the increasing depletion 
rate with increasing density). However, the time- 
averaged O2 abundance is not as strongly depen- 
dent on the shock timescale. 

5. Discussion 

5.1. Comparison with Observations 

Our results demonstrate that the time-averaged 
abundance of H2O is a sensitive function of 
the rate at which molecular gas is subjected 
to shocks. In order to properly interpret these 
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results it is useful to consider how the time- 
averaged abundances relate to observations of 
water in a giant molecular cloud (GMC). The 
gas inside a star-forming cloud can roughly be 
divided into three categories: (1) quiescent ma- 
terial that has been relatively unaffected by cur- 
rent or previous epochs of star formation; (2) 
gas that is being physically and/or chemically 
affected by current star formation; and, (3) as 
suggested here, gas that has been affected by a 
previous generation of star formation and is in 
the process of chemically and dynamically evolv- 



ing back to a quiescent state. Observations of 



H2O toward molecular material associated with 
each of these categories should find abundances 
varying from as low as x(H20) ~ 10 -6 in (1), 
and ranging upwards to < 10 -4 in (2) and (3). 
A single pointed detection of H2O emission in 
a GMC can therefore be represented by a single 
stage of the three-stage or Monte- Carlo models 
(quiescent, shock, post-shock). Because star for- 
mation will be spread throughout a cloud or core, 
the time- averaged abundance therefore refers to 
water abundances averaged over an area that 
contains both quiescent material and any asso- 
ciated gas currently being affected by local star 
formation. Thus the average abundance over a 
cloud (cloud- average) is comparable to the time- 
averaged abundances in the Monte- Carlo cloud 
model. This average abundance might not nec- 
essarily apply to an entire GMC complex (e.g. 
Orion, Gem OBI), which can extend for sev- 
eral square degrees on the sky, but may apply 
to several dense cores within a single cloud. Re- 
gions with high star- formation rates, such as 
OMC-1 which is associated with the Trapez- 
ium cluster (see Section 3), might be expected 
to have a higher probability for strong shocks, 
and therefore a high cloud-averaged water abun- 
dance. Cores with lower star-formation rates, 
such as L1641 also in Orion, would have low 
cloud- averaged water abundances. 

To test our model's ability to constrain the his- 
tory of molecular clouds it is important to ob- 



tain maps of the water emission over large spa- 
tial scales in GMC cores. Because of the strong 
absorption by atmospheric water vapor, observa- 
tions of H2O in the ISM are extremely difficult. 
As a result, the small number of detections of wa- 
ter in molecular clouds have typically been taken 
towards single lines of sight mostly containing 
luminous protostars. The very convincing detec- 
tions of water by ISO provide the greatest evi- 
dence for high water abundances in hot gas, with 
x(H20) ~ 1 — 6 x 1CP 5 inferred towards hot stars 
( [Hclmich et al. 1996| ; Ivan Dishoeck & Helmicbl 
"19961 ) an d HH 54 ( [Liseau et al. 199q ). Towards 



Sgr B2 Zmuidzinas et al. (1994) observed H^ 8 
in absorption with the Kuiper Airbourne Obser- 
vatory (KAO), while Neufeld et al. (1997) ob- 
served H2O in emission using ISO. Combining 
these observations with previous ground based 
detections Neufeld et al. (1997) estimate an 
abundance of x(H20) = 3.3 x 10 for the cooler 
outer parts of Sgr B2 and x(H20) ~ 5 x 10 -6 in 
the hot core. These observations are important 
not only because H2O was detected in dense gas, 
but also because there appears to be a range of 
water abundances. However, it is difficult to dis- 
cern whether the enhanced water abundances are 
the result of high-temperature chemistry that oc- 
curs behind shocks, high-temperature chemistry 
appropriate to gas near young stars (c.f. Doty 
& Neufeld 1997), or evaporation of grain mantle 
species (see discussion in Section 5.2) and there 
is little information on the spatial distribution of 
the water emission. 

There have been a few attempts to map the ex- 
tended emission of water. Gensheimer, Mauers- 
berger, & Wilson (1996) mapped emission of the 
quasi-thermal 3i3 — > 220 transition of H^O in 
both the Orion Hot Core and Sgr B2, and found 
that the emission originated in a compact re- 
gion (< 10"). Cernicharo et al. (1994) and 
Gonzalez- Alfonso et al. (1995) find evidence for 
widespread water emission in Orion and W49N 
using the 3i3 — > 220 masing transition of H2O. 
In Orion they argue that the water abundance is 
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quite high, :r(H20) > 10 -5 , over an extended 
50" x 50" region centered on BN-KL and the 
Orion Nebular Cluster. From these results, the 
cloud-averaged water abundance for the central 
regions of the Orion core near the Orion Nebu- 
lar Cluster is ^10 -5 which, using Figures 11 and 



12 (n# 2 = 10 6 cm 3 ; [Bergin, Snell, fc Goldsmith 



19(Q), is consistent with t s ~ 10 6 years, as sug- 
gested in Section 3. However, the determination 
of abundances from masing transitions is a dif- 
ficult task and these results must be viewed as 
suggestive until they are confirmed by mapping 
data obtained in other H2O lines. 

For molecular oxygen, which also suffers from 
strong absorption due to atmospheric O2, the sit- 
uation is even more uncertain. There exists only 
one tentative detection of 16 18 towards L134N 
by Pagani, Langer, fc Castets 1993] , which im- 
plies x(02) ~ 4 — 8 x 10~ 5 . However, a search for 



O i8 emission in different positions in L134N 



and other galactic sources by Marechal et al 
19(Q did not confirm this detection and provides 
only upper limits of02/CO < 0.1, which, assum- 
ing x(CO) = 2.7 x 10~ 4 (Lacy et al 1994), gives 
x(02) < 3xl0~ 5 . Searches, with similar negative 
results, have also been performed for molecular 
oxygen in extra-galactic sources with favorable 
redshifts (Goldsmith & Young 1989; Combes et 
al. 1991; Liszt 1992) finding x(0 2 ) < 10~ 6 . Are- 
cent study by Combes, Wiklind, h Nakai (1997) 
towards a z = 0.685 object provides the lowest 
limit to date of 2 /CO < 2 x 10~ 3 . For our 
most realistic case, the Monte Carlo gas-grain 
model, we find the time-averaged O2 abundance 
is ~ 2 x 10 -5 . Thus, the combined effects of 
shocks and gas-grain chemistry could lower O2 
abundances below the observed limits in Galac- 
tic sources, but another explanation is required 
to account for extra-galactic observations. The 
low time-averaged O2 abundance suggests that 
these results could have some bearing on the 
high abundances of neutral carbon relative to CO 
(C/CO ~ 0.1) that are observed towards a va- 
riety of star forming regions (c.f. Plume 1995, 



Schilke et al 1996). However, because the num- 
ber of dissociative shocks which destroy CO is 
not high enough, the time-averaged carbon abun- 
dance in the models is well below the observed 
value. 

The best opportunity to test these results 
will be offered by two spaceborne observato- 
ries, The Submillimeter Wave Astronomy Satel- 
l ite (SWAS) QMelnick et al. 1995Q and ODIN 
( Hjalmarson 1995| ), both of which should launch 
within the next year. SWAS and ODIN are ca- 
pable of observing and mapping the fundamental 
transition lio — > loi of H2O at 557 GHz and the 
3,3 -> 1,2 transition of 2 at 487 GHz. These 
transitions have low upper state energies (~ 26 
K) and should be readily excited in hot gas near 
star-forming sites and, more importantly, in the 
colder more extended material such as the ridge 
of dense gas in Orion (c.f. Ungerechts et al 



1997) 



Besides the computation of averaged abun- 
dances from mapping observations, combined ob- 
servations of H2O and O2 should be a powerful 
tool in examining the shock history of gas. This 
is demonstrated in Figure 16 where the water 
abundance is plotted versus the molecular oxy- 
gen abundance for the Monte-Carlo model, in- 
cluding grain depletion and desorption (Section 
4.4.2). In this plot, the majority of points trace 
an area of roughly constant water abundance of 
x(H20) ~ 10 -6 , with the O2 abundance rang- 
ing between ^(02) = 3 — 30 x 10~ 6 . This area 
has, by far, the largest number of points and 
defines the "main-sequence" of quiescent chem- 
istry. From this main sequence a shock will trace 
a line extending almost horizontally to the right, 
with constant O2 abundance. Very strong shocks 
continue at the end of the horizontal line and 
extend down almost vertically, lowering the O2 
abundance for constant H2O abundance. This 
dependence of horizontal and vertical lines oc- 
curs because the H2O is created more efficiently 
than O2 destruction (see Figures 3 and 4). This 
plot is a different way of examining the mod- 
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Fig. 16. — Points denning the plane of O2 and H2O abundances for the gas-grain Monte-Carlo model 
with t s = 10 6 years and nn 2 = 10 5 cm -3 . The spread of plotted points has been artificially increased by 
smearing them by 0.2 dex so as to give a better sense of the density of points in different regions of the 
plot. The cross represents the mean (time-averaged) abundances. 



els, as opposed to average abundances, because 
it presents the evolution in a continuous fash- 
ion; the various possible solutions define a plane 
that traces evolutionary tracks including quies- 
cent chemistry through a broad range of shock 
strengths. 

5.2. Post-Shock Chemistry and Water Ice 
Mantles 

For comparison with previous modeling ef- 
forts, our three-stage model is similar to the 
episodic models presented in Charnley et al. 
(1988a, b). However, our model is simpler in that 
we do not model the formation process of a dense 
clump nor do we account for any mixing between 
shocked and non-shocked layers. Like their mod- 
els, we find that the chemistry converges to well 
defined abundance values even when numerous 
cycles are repeated Cha.ru ley et al (1988a.) also 



fin d that the water abundance varies greatly bo 
tween shock and quiescent cycles, but they do not 
examine in detail the chemistry in the post-shock 
layer. 



An interesting result in the gas-grain models 
is that the abundance of water on grain surfaces 
in post-shock gas is quite large x(H20) gr ~ 10 -4 . 
This abundance is quite close to that inferred for 
water ice along lines of sight towards background 
stars in Taurus, where x(H20) gr ~ 8 x 10~ 5 
W'liil lei 1993| ). Thus, these results offer an al- 
ternative explanation for the large abundance 
of water in grain mantles - one that requires 
no grain surface chemistry! We stress that this 
mechanism may not be responsible for all H2O 
observed on grains. Grain surface chemistry for- 
mation of H2O must be considered - especially 
during cloud formation stages when the abun- 
dance of atomic hydrogen and oxygen is higher. 
There are also other potential desorption mech- 
anisms that could be active and are not included 
in these models, such as grain mantle explosions 
induced by UV radiation ( [Schutte fc Greenberg 



1991) or water desorption via the infrared radia- 



tion field ( Williams, Hartquist, & Whittet 1992] ). 
However, since most molecular clouds are very 
active star-forming sites, we suggest that cau- 
tion must be applied when interpreting observa- 
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tions of ice mantles in molecular clouds as the 



sol p result of grain surface chemistry. The in- 



outflows (see also Bachiller 1996; Bachiller & 
Perez Gutierrez 1997 ). These species are in- 



terpretation of high gas-phase water abundances 
towards hot stars as the result of either high- 
temperature chemistry or grain mantle evapora- 
tion of H2O is further obscured because water 
on grains may have been produced in an earlier 
shock episode. It is possible that the HDO/H2O 
ratio could discriminate between water mantles 
created in shocks or in grain mantles and we 



eluded in the three-stage chemical model pre- 
sented in Section 4.3, but we do not present the 
results in detail here because of uncertainties in 
the high-temperature reaction rates. However, 
these results do have some bearing on the chem- 
istry of these species in the post-shock layer. The 
formation pathway of CH3OH in the gas phase 



is linked to a reaction of CH3 with H2O (Millar 



ar ej in the process of investigating this question Herbst, fc Charnley 199l| ) . Thus, it is possible 

(Bergin, Neufeld, & Melnick 1997). 

The large abundance of water on grains, with 
CO remaining in the gas phase, will also al- 
ter the ratio of total carbon-to-oxygen (C/O) 
in the gas phase. Comparison of chemical the- 
ory with observations of molecular abundances in 
GMC cores associated with massive star forma- 
tion has found that the chemical abundances of 
many species are best reproduced with C/O ra- 
tios greater than the solar value (> 0.4; |Blake et 



that the enhanced abundance of water through 
high-temperature reactions could lead to larger 
CH3OH abundances through this reaction. In 
the three-stage model, the CH3OH abundance in 
the post-shock layer is x(CH30H) = 1 — 5 x 10~ 8 
when the shock temperature ranges from 1000 to 
2000 K. These abundances are at least an order 
of magnitude below values inferred in molecular 
outflows, which can be as high as x(CH30H) ~ 
10~ 6 ( Bachiller et al. 1995 ). The mechanism of 



al. 1987; Bergin et al. 199*7] ). For chemical mod- 
eling this is the result of depletion of the initial 
abundance of oxygen relative to carbon, which 
reduces the abundances of small oxygen-bearing 
species that are major carbon destroyers in the 
gas phase. In the three-stage gas-grain model 
shown in Figure 7, the C/O ratio in the post- 
shock gas at t ~ 10 5 years is C/O ~ 0.7. It is dif- 
ficult to gauge the effect of this on other species, 
because of the lack of certainty with regard to 
high-temperature reactions. It is therefore pos- 
sible that the C/O ratios inferred in Blake et al. 
(1987) and Bergin et al. (1997) are indicative of 
the core formation process, which could involve 
gas temperatures rising high enough to convert 
atomic oxygen to water. In this case, these re- 
sults suggest a simple mechanism to place large 
amounts of oxygen on grain surfaces while still 
keeping most of the carbon in the gas-phase. 
Other molecular species, such as SO, SiO 



flMartin-Pintado, Bachiller, fe Fuente 1992[) , and 
CH3OH flBachiller et al. 1995Q have been ob- 
served with enhanced abundances in energetic 



methanol enhancements may therefore be the re- 
sult of grain mantle evaporation or an uniden- 
tified high- or low-temperature pathway. The 
abundance of CH3OH accreted onto grain sur- 
faces provides another constraint. In our model 
the abundance of CH3OH ice in the grain mantle 
is only 0.1 percent of the water-ice abundance. 
This ratio is below that observed towards NGC 
7538 or W3 3A where x(CH 3 0H)/x(H 2 0) ~ 10- 
40 percent ( Allamandola et al. 1992] ). Our mod- 
els do not set constraints on the chemistry of SO 
and SiO because the large abundances of these 
species in outflows may be the result of sputter- 
ing of grain refractory and/or mantle material 
(c.f. Schilke et al. 1997; Caselli et al. 1997). 

When the abundances of these species are en- 
hanced, through any mechanism, the high abun- 
dances will persist until the timescale for the 
individual molecule to deplete onto grain sur- 
faces. If the dust temperature is high enough 
to keep a given species in the gas-phase, or for 
pure gas-phase chemistry, the lifetime will be 
T ra = 4 — 7 x 10 5 years. In these models all 
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SiO depletes onto grain surfaces at Td ep ~ 10 4 
years (for n# 2 = 10 5 cm -3 and assuming a stick- 
ing coefficient of unity). For CH3OH abun- 
dance enhancements we find that the depletion 
timescale for gas-grain chemistry is equal to that 
of H 2 (T dep (H 2 0) ~ 10 5 years). The dispar- 
ity in post-shock lifetimes between CH3OH, H2O, 
and SiO suggests that differences might exist be- 
tween younger outflows, perhaps those associ- 
ated with so-called Class sources, and the more 
evolved sources (e.g. Class I). A survey of such 
sources in these tracers might prove to be useful 
in probing the links between evolutionary effects 
observed in outflows and the relationship to the 
driving source. Another possibility is that differ- 
ences could also be found within different com- 



ponents inside a single outflow (see Bachiller & 



Pe rez Gutierrez 1997 ). 

The abundance of some molecular species, no- 
tably HCO + and N2H + , are adversely affected 
by the high water abundance. In Figure 7, for 
the three-stage model, the abundance of HCO + 
is decreased through reaction R6 when the water 
abundance is raised. Thus, the abundances of 
these two important molecular species should be 
anticorrelated. Because N 2 H + also reacts with 
H2O, a similar anticorrelated behavior would be 
found between the abundances of N2H + and H2O 
and may account for the low HCO + and N2H + 
abundances in hot regions such as the Orion Hot 
Core (c.f. Blake et al. 1987) and in the LI 157 
outflow (|Bachiller Perez Gutierrez 1997| ) . 

6. Summary 

We have investigated the use of simple oxygen- 
bearing molecules, H2O and 2 , as tracers of the 
shock history of molecular gas. To this end, we 
have constructed a coupled chemical and dynam- 
ical model of dense molecular gas that is sub- 
ject to heating by non-dissociative C-type shocks, 
and subsequently cools. To constrain the effects 
of MHD shocks on molecular gas we use the shock 
models of Kaufman & Neufeld (1996a,b). We 
have constructed two models and concentrate, in 



particular, on the chemistry in a layer of post- 
shock cooled gas. In the first model we exam- 
ine the chemical evolution in three-stages: (1) as 
molecules form and increase in abundance within 
quiescent, pre-shock gas; (2) the rapid chemical 
changes precipitated by the passage of a C-shock; 
and, (3) the post-shock chemical evolution that 
occurs as the gas cools. The results of the three- 
stage model are combined with observational es- 
timates of shock rates within molecular gas in re- 
gions of molecular outflows. A second, more gen- 
eral, model uses a Monte-Carlo method to exam- 
ine the chemistry in a parcel of gas that is subject 
to random shock heating over a cloud lifetime. 
Because water is also an abundant molecule on 
the surfaces of interstellar grains, we have sep- 
arately examined the pure gas-phase chemistry 
and the gas-phase chemistry including grain sur- 
face molecular depletion and desorption. 
The principal results from this study are: 

1) For v s > 10 km s _1 , high temperature 
neutral-neutral chemical reactions are efficient at 
converting all of the available oxygen into H2O 
within the post-shock region before the gas has 
an opportunity to cool. For lower shock veloc- 
ities, the water abundance will not change sig- 
nificantly from its quiescent values due to the 
rapid post-shock cooling to temperatures less 
than those required to initiate these neutral- 
neutral reactions. This result is in agreement 
with previous investigations. Similarly, for v s > 
26 km s _1 , all O2 will be destroyed inside the 
shock within the cooling time. 

2) For pure gas-phase models, the enhanced 
H2O abundance behind a shock will decrease 
and the low-temperature chemistry will be re- 
asserted at r ra ~ 4 — 7 x 10 5 years. As a result 
of cosmic-ray driven chemistry, r ra is nearly in- 
dependent of the gas density. For models which 
include gas-grain interactions, a density depen- 
dence is found as the large abundance of water 
in the gas-phase will decrease though accretion 
onto grain surfaces at the depletion timescale. 
The depletion timescale has an inverse depen- 
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dence on the density and is less than r ra for n# 2 
> 10 5 cm~ 3 . 

3) We find that the time-averaged abundance 
of H2O is a sensitive function of the shock fre- 
quency. This abundance is approximately the 
weighted average of the amount of time gas 
spends in the pre-shock, shocked, and post-shock 
stages respectively, and determines the average 
abundance of water mapped over a large region. 
Thus these models predict that the abundance of 
H2O, and to a lesser extent O2, can be used to 
constrain the physical history of molecular gas. 
The observations required to test these predic- 
tions can be performed from The Submillimeter 
Wave Astronomy Satellite and ODIN. 

4) The abundance of water-ice on grain sur- 
faces can be quite large in the post-shock gas 
layer. The abundance is often x(H20) gr > 10~ 4 
and is comparable to that observed in molec- 
ular clouds. This offers an alternative method 
to create water-ice mantles without resorting to 
grain surface chemistry - passage of a C-type 
shock, followed by rapid hydrogenation of gas- 
phase oxygen, and subsequent depletion of water 
onto cold dust grains. Despite the large amount 
of water trapped on grains, the gas-phase chem- 
istry is maintained because CO and O2 do not 
deplete. This differential depletion of H2O com- 
pared to CO and O2 can also effectively raise the 
C/O ratio in the gas phase to values that are in- 
ferred from previous observational and theoreti- 
cal comparisons of chemical abundances in GMC 
cores. 

5) The abundance of a few molecular species 
are adversely affected by the high water abun- 
dances in shocked regions. In particular the 
abundances of N2H + and HCO + are predicted 
to have abundances that are anticorrelated with 
changes in the water abundances. 
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